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Earlier work on dynamical critical phenomena in the context of magnetic hysteresis for uniaxial 
(scalar) spins, is extended to the case of a multicomponent (vector) field. From symmetry arguments 
and a perturbative renormalization group approach (in the path integral formalism), it is found that 
the generic behavior at long time and length scales is described by the scalar fixed point (reached for 
a given value of the magnetic field and of the quenched disorder), with the corresponding Ising-like 
. exponents. By tuning an additional parameter, however, a fully rotationally invariant fixed point 

can be reached, at which all components become critical simultaneously, with O (n)-like exponents. 
Furthermore, the possibility of a spontaneous non-equilibrium transverse ordering, controlled by 
a distinct fixed point, is unveiled and the associated exponents calculated. In addition to these 
central results, a didactic "derivation" of the equations of motion for the spin field are given, the 
scalar model is revisited and treated in a more direct fashion, and some issues pertaining to time 
dependences and the problem of multiple solutions within the path integral formalism are clarified. 
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I. INTRODUCTION 



In a great variety of non-equilibrium situations, critical behavior is observed as a system evolves from one of its 
possible states to another. Some examples are charge density waves, fluctuating interfaces and lines, cracks and 
fractures, and the Barkhausen effect in magnets. These systems evolve, respectively, from a state without current to 
a state with current, from a stationary to a moving state, from a connected to a ruptured state, from a downward 
to an upward magnetization. Under specific conditions (i.e. preparation of the system), the transition between 
the two states is critical (or continuous), exhibiting diverging correlation lengths, and scaling laws. The qualitative 
descriptions of the dynamics of the different physical situations mentioned are very similar in the key parameters and 
mechanisms which govern criticality. The quantitative descriptions are also close, in that the dynamics of motion can 
be described by continuum field equations, and share many common features. In this article we focus on magnetic 
systems. More specifically, on a lattice of spins with ferromagnetic exchange (coupling). While our qualitative and 
quantitative analyses will be in this framework, some aspects of the discussion may apply to other systems, such as 
£ — . ■ depinning transition of flux lines or fractures in disordered media. 

What do we mean by "the key parameters and mechanisms which govern criticality"? Consider a ferromagnet in 
an external magnetic field which increases slowly from — oo to +oo. Each spin feels a local field equal to the average 
of the surrounding spins multiplied by the coupling constant (JMi in the case of the i th spin), plus the external 
magnetic field H . At H — — oo, all the spins point "downward" (Mj = —1 for unit spins), thus JMi = — J initially. 
At zero temperature, each spin simply points in the direction of the local field JMi + H, and so none of the spins 
change before H reaches J, at which point they all flip upwards. The magnetization M thus jumps from —1 to +1 
rj^2 ' &t H — J . This scenario for a perfectly clean system is modified by introducing some disorder. At each lattice point 
occupied by a spin, add a random field, hi, to the local field, JMi + H. (Specifically, let hi be an uncorrelated random 
variable, chosen from a Gaussian distribution centered at zero.) Then, the spins flip in a much less coherent way: as 
soon as JMi + H + hi becomes positive, the i th spin flips. The upward h^s enhance the increase in magnetization 
^ , for low H, whereas the downward his suppress it for high H. This results in the reduction of the magnitude of the 
jump in magnetization. Clearly, if we broaden the random field distribution, i.e. increase the amount of disorder, 
the discontinuity in M is further suppressed, until the curve M(H) eventually becomes smooth for a high enough 
disorder (Fig. 1). We can imagine a sequence of hysteresis curves, corresponding to a succession of increasing amounts 
of disorder, say for the variance of the random field, h 2 , going from zero to infinity. The curve displays a discontinuity 
for small h 2 , and is smooth for large h 2 . The transition between the two regimes occurs at a critical point reminiscent 
of continuous or second-order phase transitions : at the critical amount of disorder the discontinuity collapses to 
a point, at which the slope is infinite. This is referred to as a critical hysteresis, for which, at a given magnetic field, 
the susceptibility diverges. The amount of disorder and the magnetic field are the two parameters we have to tune to 
observe criticality. 
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FIG. 1. Schematic hysteresis curves for different values of the disorder 7?. Left: R < R c (discontinuous hysteresis); center: 
R — R c (critical hysteresis); right: R > R c (smooth hysteresis). In each case, the lower (upper) curve corresponds to an 
increasing (decreasing) magnetic field. 



In the last few years, disorder-induced critical hysteresis in magnets has been the subject of much interest [gj — jSJ] - 
Dahmen, Sethna, and others studied this problem via a mean-field approximation, one-loop momentum-space renor- 
malization, and numerical simulations. Also, they describe a mapping of this non-equilibrium problem onto the 
equilibrium random field Ising model, which can in turn be mapped (close to the upper critical dimension) onto the 
pure Ising model in two lower dimensions |t]]. Throughout their work, they consider a scalar order parameter. They 
study the dynamics of an Ising (or discrete) spin field driven by an increasing magnetic field and in the presence of a 
random field, at zero temperature. 

The question we ask here is the following: How are the phase diagram and critical behavior modified if the order 
parameter is vectorial instead of scalar? More precisely stated: in the renormalization group (RG) framework, is 
the fixed point which controls the above-mentioned hysteretic criticality, one and the same for both the Ising and 
the vectorial cases? And if not, how do the exponents differ? In equilibrium, the disordering of scalar and vector 
systems is described by distinct universality classes @]. Also, in the closely related context of depinning transitions, 
the distinction between interfaces (scalar) and flux lines (vector) was noted in Ref. || . 

The answer can readily be guessed on symmetry grounds. Indeed, symmetry considerations lead to two distinct 
cases. In the first, and generic one, the critical hysteresis curve is such that the susceptibility diverges at a non- 
vanishing value of either the magnetic field or the magnetization. Then, although we consider continuous spins, the 
full rotational symmetry of the problem is broken at the critical point, a unique preferred direction is picked, and 
Ising-like critical behavior results. It is similarly argued, in Ref. pi, appendix E, and Ref. S, appendix L, that the 
universality class of the random field scalar model extends also to random bonds scalar models (with a positive non- 
zero mean of the bonds' values) and to random anisotropics O (n) models. On the other hand, if both the magnetic 
field and the magnetization vanish when the susceptibility diverges, we may have a fully rotation invariant system 
at the critical point. In that case, we expect O (n)-like criticality, with exponents that differ from those of the Ising 
model. Evidently, a vanishing magnetization at H = is not a sufficient condition for a full rotational invariance. In 
particular, due to its history, the system might well display higher order anisotropics, such as, e.g. s| ^ s j_ Q ? where 
sii is the component of the spin field parallel to H and sj_ a ia any parpcndicular component. This issue is resolved 
by a renormalization group analysis, which confirms our various guesses. Furthermore, it discloses the possibility of 
a "transverse critical point" , corresponding to an instability of the magnetization component perpendicular to the 
external magnetic field. 

The present paper is organized as follows. In Sec. II, we construct the equation of motion of a vector spin field. 
A path integral formalism is described in Sec. Ill, with which the problem of renormalizing the equation of motion 
is recast into that of renormalizing a partition function. The time dependences, as well as the subtleties associated 
with "many energy minima" , are also examined. In Sec. IV, the renormalization group treatment of the problem is 
presented. First, we define the coordinates' and fields' rescalings, and calculate the free propagator. Then, we discuss 
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successively the scalar and vector models. For the latter, the different cases (hysteretic or non-hysteretic, longitudinal 
or transverse criticality) are analysed, and the corresponding recursion relations and exponents are obtained. 



II. EQUATIONS OF MOTION 

The equations of motion for a scalar field are discussed in Refs. and their generalization to a multicomponent 

field is straightforward. Nonetheless, for completeness and to emphasize our perspective, we present here a didactic 
introduction to the equations of motion for vectorial spins, at zero temperature Consider a (/-dimensional lattice, 
with a spin G 5f n at each site i, subject to a magnetic field which changes slowly from — oo to +00, say linearly 
in time, H = fit. The rate fi can be made arbitrarily small in magnitude and points along the first axis of our 
coordinates, i.e. Hi = Qt = H, H2 = ■ ■ ■ = H n = 0. The time dependent magnetic field implies a time dependent 
energy function TL (t). At zero temperature (T = 0), the spins simply follow the local minimum of this energy function 
according to 

■qd t s l = - — , (1) 
OS, 

starting from a uniform downward pointing configuration at t — —00. The parameter rj controls the relaxation rate 
of spins towards the time-dependent local energy minimum of TL. The smaller 77, the faster spins relax, and the less 
they lag behind the energy minimum poj . 

The following "key features" guide us in constructing the Hamiltonian TL. First, to describe a ferromagnet, we 
include in TL couplings Jy which tend to align the spins. In addition to the external uniform field H which drives the 
system, we include quenched random fields hj. The hj's are uncorrelated Gaussian random variables, chosen from 
the distribution 



p[h]=iVexp -V-i , (2) 




where A is a normalization factor. For calculational convenience we shall work with soft spins (whose magnitude can 
take any real value), which can be thought of as a coarse-grained picture of a hard spin field. To avoid the unphysical 
instability of spins diverging in magnitude, we introduce an on-site potential V (sj), which constrains the magnitude 
to remain close to 1 (or some finite number). The potential V is spherically symmetric (a "double well" in the scalar 
model and a "Mexican hat" in the vector model) and is expressed through its Taylor expansion about the origin, as 

V = s?-|( Sl 2 ) 2 + .... (3) 

Whether or not V is analytic at the origin is unimportant, since |sj| is constrained to be close to 1 (and not to 0). 
The full Hamiltonian is now given by 

TL = -- J] JijSi [-H • s, - hi • Si + V (a,-)] . (4) 

i,3 i 

The gradient descent with this Hamiltonian leads to the equation of motion 

■qdtSj = JjjSj + H + hj - — . (5) 

j 

Assuming that Jjj is a function of the separation between spins results, in the continuum limit, in 

r]d t s (x) = J d d x'J (x - x') s (x') + H (t) +h (x) (6) 

While rewriting the problem in the continuum limit, we must impose some limit on how fine-grained the spin field 
s (x) may be because of its lattice origin. In other words, s (x) is a superposition of Fourier components whose wave 
numbers are restricted from zero to some cutoff A. 

Finally, if the exchange function decays fast enough (as its argument increases) for its Fourier component to be 
non-singular at the origin of momentum space, i.e. if J (q) = 1 — Kq 2 + • • •, then 
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J d d x'J (x - x') s (x') « s (x) + ^V 2 s (x), (7) 
and the equation of motion can be written as 

8V 

nd t s = KV 2 s + H + h-—, (8) 

as 

where the coefficient c\ in the the expansion of the potential has been modified in order to take the s (x) term of 
Eq. (Q) into account. 

III. PATH INTEGRAL FORMALISM 

In order to identify the critical properties of our model, we should ideally solve its equation of motion. In practice, 
we study the behavior of Eq. (||) under a coarse-graining transformation. This allows us to locate and characterize a 
scale-invariant point. As a first step, we recast the equation of motion into a path integral (or generating functional), 
which incorporates the whole history of the system. The generating functional is written as the sum over all paths 
of the exponential of some action, which is then renormalized perturbatively. The advantage of reformulating the 
problem in this way is that we can express the perturbative treatment in a diagrammatic fashion similar to other field 
theories. 

We define the generating functional pJ| simply as the sum over all paths of a delta function which makes each spin 
follow the time evolution given by Eq. (pft, i.e. 

Z = j Vs 8 (s - solution of Eq. (8)) 

= J Vs 5 ^-r)d t s + KV 2 s + H + h - x (Jacobian) , (9) 

where Vs stands for Yl {over "all t" , "all x" , a = 1, • • • , n} ds a (x, t). The Jacobian merely normalizes the value of Z 
to unity, and we shall henceforth ignore it jl2] . Let us rewrite the delta function in its representation as the integral 
of an exponential [2ir6 (/) = J (is exp (isf)]. Then, after absorbing a factor i in a redefinition of s, and dropping an 
infinite multiplicative constant (along with the Jacobian) , we have 

Z = J Vs Vs exp I J dtd d x s • ^-J?d t s + K V 2 s + H + h j 

ee J VsVs exp (S) . (10) 

The generating functional Z enables us to evaluate all correlation and response functions. For example, the solution 
of Eq. (1) is 



s sol 



(x, t) = J Vs Vs s (x, t) exp (S) , (11) 



and its response to the magnetic field is given by 

SsT 1 (x,i) 



JvsVsJ d d x' si (x', t') si (x, t) exp (S) . (12) 



SH (*') 

Also, we can change the origin of time by a trivial reparametrization of the magnetic field, as e.g. in 

s so1 (x, t + e; H (t)) = s so1 (x, t; H (t + e)) , (13) 
The latter expression takes the form s so1 (x, t; H (t) + fie) if H is increased linearly in time at a rate f2, whence 



s- (x, t + c) - s« (x, t ) . / VsVa . (x, t , e x p ( S ) {«p (/ *W h0t ) - 1} , (u, 
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from which it follows that the dynamic susceptibility is calculated as 



ds so1 (x, t) 
dt 



= Q J dt'd d x' J Vs Vs s x (x', t') s (x, i) exp (S) . 



(15) 



Since we are interested in the average of the correlations and responses over the random field, from now on we deal 
with the average of Z. This enables us to forget the stochastic variable h a , trading it for a new term in the "averaged 
action" . Taking advantage of the Gaussian nature of h a , 



Z = J VsVs ex P 1^ dtddx § • \ + KV 2 s + H - j j exp J dtd d x s • h 

= J VsVs exp (S) , (16) 

with (using Eq. (g)) 

S = J dtd d x s • (^-r)d t s + KV 2 s + H ^ j + j dtdt'd d x (x, t) ■ s (x, t') . (17) 

We have reformulated the theory, originally described by a dynamical differential equation, in terms of an action 
5 [s (x, t)], which depends on the entire history of all spins. Thus S is a functional of the path which the system 
follows; the probability weight exp (5) picks the physical path and averages it over disorder. We can then study the 
symmetries and renormalization of the theory, as for equilibrium field theories. 

The motion comprises two time scales: r\ and (dH j di)~ x . The behavior of the system depends of course on the 
ratio of the two, and not on their respective values. For the calculation of static exponents, we let 77 — s- 0. Consider 
for example the exponent with which the correlation length diverges. A diverging correlation length gives rise to 
an infinite susceptibility, detected by a non-vanishing response (of the magnetization) to an infinitesimal increase of 
the magnetic field. Clearly, such a behavior is obtained in our problem only if the magnetic field increases infinitely 
slowly [i.e. (dH/dt) —*•()], or equivalently, if 77 ^ 0. As 77 is a measure of how much the system lags behind its local 
energy minimum, for rj — > the system always stays at the time dependent minimum (as seen by setting 77 = in the 
equation of motion). In other words, the spin avalanches resulting from a small change in H spread instantaneously 

In the 77 — > limit, time evolution of the spin field is simply motion with the minimum in the energy landscape. 
There is nevertheless a subtlety involved as to the choice of the minimum. For example, because of the double-well 
shape of the potential, a scalar spin has the choice, during some time interval of its history, between two positions, 
both of which locally minimize the Hamiltonian. If we solve the equation of motion, there is no ambiguity, since we 
specify one of the energy minima as the initial condition, and then follow its time evolution. In particular, the initial 
condition corresponding to the case of a magnetic field which increases from -co to +00 is one in which all the spins 
occupy the left minimum of their double wells fl4f| . In the path integral formalism, however, no initial condition is 
specified. The weight e picks all possible solutions, corresponding to different initial conditions. Let us illustrate this 
point with, as in Ref. JL5L the zero-dimensional non-random version of our model, defined by the equation of motion 

r)d t s = -^-(-Hs + as 2 + bs 4 ) , (18) 

with 77 — > 0. We can imagine s as a bead sitting at the minimum of the quartic potential Q (s) = —Hs + as 2 + bs , 
and moving with it. With a < and b > 0, Q (s) has a single minimum if \H\ is larger than some value A, and 
two minima if \H\ is smaller than A. For H very negative, the bead sits at the single minimum of the curve, which 
becomes the left minimum when H is between —A and A. At H = A, the left minimum disappears and the bead 
moves to the right minimum. Thus, in particular, s (-H) ^ ~s (H) and s (H — 0) 7^ 0: the motion of s is hysteretic 
(Fig. 2). 
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FIG. 2. Metastable solutions of the single-spin system. The physical solution for an increasing magnetic field (solid line) 
follows the lower curve (s_) as long as it is present, then jumps along the dashed line to the upper curve (s+). 



The action for this model is 

S = J dts (-rjd t s + H-as- 6s 3 ) . (19) 

Let (s) be the average of s with respect to the weight e . Calculating (s) perturbatively, it is easily seen (by counting 
the possible occurrences of s and s in a diagram) that each diagram comes with an odd number of 77's. Therefore, 
(s (-H)) = — (s (H)), and in particular, (s (H — 0)) = 0, in apparent contradiction with the above solution. But 
as mentioned earlier, the partition function is merely the integral of a delta function which imposes the equation of 
motion. In the present case, it imposes 

d 

— (-Hs + as 2 + 6s 4 ) = 0. (20) 

Thus, in calculating (s), we pick all minima of the quartic form. In other words, (s) is the sum of two terms, 
(s) = (s_) + (s+), corresponding to the left and right minima. The physical, "bead" solution is equal to (s_) up to 
H = A and then (s + ) for H larger than A. Similarly, in the case of our original problem, the quantity s so1 (x, t) of 
Eqs. ( pT| ) and ( |l2| ) does not coincide with the physical solution we are looking for. In addition to the latter, s so1 (x, t) 
contains other unwanted solutions corresponding to additional energy minima. We shall come back to this difficulty 
and circumvent it in the next section, both for the scalar and the vector models. 



IV. PERTURB ATIVE RENORMALIZATION 



A. Coordinates and Fields Rescalings, the Free Propagator 

Our renormalization group transformation consists of the usual three steps. First, we coarse-grain the system, i.e., 
integrate out the modes with wave number between A/6 and A (6 > 1). Second, we rescale coordinates (i.e. change 
our length and time units), by setting 

x -> 6x, t -» b z t, (21) 



G 



or, equivalently 



q^6 x q, u^b~ z uj. (22) 

With this change of units, the coarse-grained fields vary on the same length scales as the original ones, and the lattice 
cutoff is preserved. Third, we rescale the fields according to 



or, equivalently 



With the choice 



s (6x, b z t) -> (s (x, t) , s (6x, b z t) -> C§ (x, t) , (23) 

b (6 _1 q, 6 _ *w) ^6 d+z Cs(q,o;), § (ft^q, b- z Lu) -> 6 d+2 C§ (q, w) . (24) 

2 = 2, (25) 

C = 6 2 -*, (26) 

C=&- 2 "*, (27) 

the time derivative and Laplacian terms of the action, as well as the ss terms, become scale invariant. The recursion 

relations will be calculated to the lowest non-trivial order in the interaction. To this order, no corrections of the 



parameters 77, K, or R of the action in Eq. (17) occur in the coarse-graining transformation, i.e. 77, K, and R arc 
invariant under our perturbative renormalization. 

As in the momentum space RG treatment of the </> 4 theory, we consider the quadratic part of the action as a Gaussian 
(free) theory, and the rest as a perturbation (interaction). Following Refs. PJ-0], for calculational convenience we treat 
the disorder-induced ss term as an interaction, instead of including it in the Gaussian part. The free theory thus 
consists only of the ss part of the action, and the corresponding bare propagators are 

(s a (q,u)sp(q',u}')) = 

(s a (q,uj)sp(q',Lu')) = , (28) 

(Mq,^) (q>')> = S a p2nS [u> + lo') (2irf S d (q + q') _^ Kqi _ ra 

where (-) denotes an average with respect to the Gaussian weight, and the indices run from 1 to n. The parameter 
r a is the q-independent coefficient of the quadratic s a s a term in the action. Fourier transforming back in time, we 
have 

( if t' < t 

(s a (q, t) s fj (q', f )>„ = I 5af} {2n) d sd (q + q/) cxp [_^!^ y _ t) ] h if t , > t ■ (29) 

(This is calculated for r a < 0. As we shall see, r a is indeed negative at criticality.) In the 77 — > limit, the propagator 
becomes 

(S a (q, t) sp (q', O> = 5 af} (2nf 5 d (q + q') - ±- S (f - t + ) . (30) 

That is, the contraction of s a (t) and s a (t') is non-vanishing only if the two times are equal (actually, only if t' is 
infinitesimally higher than t). With this propagator, it is easily seen diagrammatically that, although the disorder- 
induced ss terms couple different times, a renormalized vertex at time t is a function only of the other vertices at 
the same time t. Thus, slices of the action at different times renormalize independently from each other, and flow 
to their respective fixed points. This justifies the procedure of Refs. [0-^) of setting H constant for the calculation 
of static exponents. (In what follows, we shall need to correct the action of Eq. (JTt]) to take care of the problem 
of multiple solutions. We shall, for example, expand S about a uniform but time-dependent value of the field, say 
some function a(t). Hence the vertices (coefficients in S), and in particular the "masses" r a , become functions of 
the parameter a (t). Note that the above analysis of the 77 — > limit is then legitimate only if a (t) is continuous, so 
that c(i + ) = cr(t). Below, we shall define our ct's in terms of the magnetization or analogous quantities. Thus, our 
analysis holds if we approach criticality from the high disorder side.) 
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B. The Scalar Model Revisited 



The zero-dimensional model discussed in Sec. Ill is nothing but the single-spin equivalent to the scalar field. The 
apparent contradiction mentioned therein is also present in the full model. Let us call s_ (x, t) and s+ (x, t) the 
solutions of Eq. (p|), for a magnetic field H increasing from -co to +00, or decreasing from +00 to —00, respectively. 
The magnetization measured experimentally is the average over space, or equivalently over the random field, of the 
above solutions, 



m±(fl-(t)) = s±(x,t). (31) 

Generically, as is inferred from the single-spin case and observed experimentally, the magnetization displays hysteresis. 
In particular, m± {—H) ^= —m± (H) and m± (H = 0) 7^ 0. On the other hand, the action of Eq. ( p"7| ) is invariant under 
the transformation (s,s,H) — ► (— s, — s, — H). This implies that the average of s so1 [Eq. ([Tl|)] satisfies s sol (—H) = 
— s so1 (H), and s so1 (H = 0) = 0. As explained above, s so1 contains unphysical contributions (corresponding to the 
many minima in the energy landscape) in addition to the physical solution (s_ in the case of an increasing H), i.e. 

s so1 (x, t) = s_ (x, t) + s' (x, t) + s" (x, t) + ■ ■ ■ . (32) 

Hence the action 

S = J dtd d x s (-r]d t s + KV 2 s + H + as + bs a + • ■ ■) + J dtdt'd d x (x, t) s (x, t') , (33) 

does not describe the magnetization, unless it is corrected in such a way as to remove the unphysical solutions. These 
are inopportunly introduced through Eq. (^), which should in fact be written as 

J Vs 5 ^-r]d t s + KV 2 s + H + h- x (Jacobian) 

Vs {5 (s - s_ (x, t)) +S(s-s' (x, f)) + 5 (s - s" (x, t)) + • • •} . (34) 



Therefore, if we substract the quantity S (s — s') + S (s — s") + ■ ■ ■ from the functional e s , we obtain a well defined 
theory, which incorporates only the physical solution and yields all the correct correlation and response functions. 
Although this method is impossible to implement (since it explicitely involves the unphysical solutions), a restricted, 
weaker version is applicable, and fully serves our purposes. The investigation of the critical behavior and the calculation 
of the corresponding exponents relies primarily on correlation and response functions, namely the field density s and 
its susceptibility ds/dt, which are linear in s. For such objects, 5 (s — s_) + S (s — s') + S (s — s") + ■ ■ ■ can be replaced 
with S (s — s_ — s' — s" — •■•). Furthermore, only averaged quantities are of interest to the study of the critical point, 
and linearity allows to perform the average inside the argument of the delta function, leading to 

S (s - m- (*) - cr (i)) , (35) 

where 



a{t) = s'(-x,t) + s"(-K,t) + --: (36) 

A comparison with Eqs. (^) and ( |l6| ) implies that a weight e s which properly describes the averaged theory is obtained 
by adding a (t) to the argument of the action. Indeed, for 

(s) = J VsVss(x,t)e s[S ' s] , (37) 

we have (s) — m_ (t) + a (t); therefore shifting the field by a(t) yields 

Vs Vs s (x, t) e s[s ' s+a] = JvsVs [s(x, t) - a (t)] e s[S ' s] = (s) - a (t) = m_ (t) . (38) 

The average spin dynamics is thus properly described by the corrected action 
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S = S [a, s + a] 

= J dtd d xs(-r]dts + KV 2 s + H + A + A lS + A 2 s 2 + A 3 s 3 + ■••)+ / dtdt'd d x (x,f) s (x,i') , (.'!!)) 



R. 



where 



' A — —r\d t o + aa + ba 3 

A 1 = a + iba 2 H 

A 2 = 36cr H 

A 3 = & + -.. 



(40) 



Although we cannot obtain the precise form of a (t) without solving for the many minima, its qualitative shape is 
easily found. Consider a very large (positive or negative) magnetic field. The unphysical minima are due only to the 
spins in a very large random field h, such that the potential they feel still has two minima. When we take the average, 
the minima s' + #" + ••■ contribute to a only with a very small weight [oc exp (-h 2 /2R)}, implying 

a -> ± , as H -> ±00. (41) 

Furthermore, it follows from the definition of a that a = — m_ at H = 0. Thus, a (H) has a bell shape, shifted 
to the right since the sum of a and is an odd function of H (Fig. 3). The parameter a (t), which corrects the 
coefficients in S, breaks the up-down symmetry present in the action of Eq. ([53|). This is physically required, since in 
an hysteretic system, the (non-equilibrium) magnetization breaks that symmetry, in particular at H = 0. 




FIG. 3. Schematic shape of the parameter a as a function of the magnetic field. 



To discuss the relevant terms in the action, we apply a renormalization transformation in 6 — e dimensions. All 
terms of the action higher than cubic in s are irrelevant, and after a large enough number of coarse-graining steps, 
we are left with the renormalized action 

S ' = J dtd d x s (-r]dts + KV 2 s + A + A x s + A 2 s 2 + A 3 s 3 ^) + J dtdt'd d x (x, t) s (x, t') , (42) 

where Aj>i's are functions of H through a, and A 3 is negative to prevent the spins from diverging in magnitude. 
Furthermore, it is easily shown that A; is an even (odd) function of a for i odd (even). In particular, since a — > as 



9 



if — ► ±00, we also have that A2 — > as if — » ±00. Now, in order to find the critical point, let us expand the field 
about some value /x (t), so as to cancel A2. With the choice /i = —^2/3^3, the action in terms of s' = s — fi becomes 

S = J dtd d x s [-r,d t s + KV 2 s' + A' + A[a' + i 3 s' 3 ) + / dtdt'd d x (x, t) s (x, t') , (43) 

where the coefficients have been corrected by fi. Since (s) — > ±00 as if — > ±00, (s) crosses /j at a given ifo, i.e. 
(s (ifo)) = /U (^o)- Hence (s' (ifo)) = 0, which implies in general that A' (ifo) = 0. The action thus reduces to that 
studied in Refs. pHl, and is critical for a specific amount of disorder R. 



C. The vector model 



1. The appropriate action 

The first question in the case of vector spins is whether the action of Eq. (|l7|) correctly describes the system. The 
single-spin problem is identical to that of a bead moving in a Mexican hat tilted by H + h, the sum of the external 
H (t) and the quenched random field h. The bead simply turns around the bump of the hat, i.e. the spin always points 
in the direction of H + h, as in the equilibrium problem. The motion is governed by a single minimum, and there is 
no hysteresis. Similarly, mean field theory reduces the time evolution of the spin field to that of a single degree of 
freedom, and yields no hysteresis for any value of the disorder. These observations may suggest that Eq. ( |l7| ) properly 
describes the problem, and that criticality occurs at H = 0, simultaneously for the components of the field parallel and 
perpendicular to H. However, although a single spin has only one minimum in its energy landscape, a configuration 
of several spins may have many. Consider a system composed of two spins of unit length, in an increasing magnetic 
field H, and with random fields +A and —A perpendicular to H, as illustrated on Fig. 4. 





H 



FIG. 4. Illustration of the two-spin toy model described in the text. 



The Hamiltonian for this system is 

H = -Jsi ■ s 2 - A- (si - s 2 ) - H • (si + s 2 ) . (44) 

If 81 and 62 are the angles of the spins with respect to H, the minimum energy path that the spins follow imposes 
81 = —82 = 8, and in terms of 8, the energy is 

H = 2 jj(sin(9) 2 - A sin 8 + if cos 6>j - J. (45) 

For A < 2J and if = 0, the energy as a function of 8 is a symmetric double well centered on 8 = ir/2 (Fig. 5), 
similar to the scalar single-spin energy landscape, with two minima at sin (9 = A/2 J. A non- vanishing ff tilts the 
double well, and ultimately suppresses one of the two minima. The ferromagnetic interaction causes the two spins 
to pull each other back before jumping ahead, thereby investing the n-component field's motion with a hysteretic, 
uniaxial-like character. That is, the two-spin toy system as a whole goes over an energy barrier, reminiscent of the 
motion of a scalar spin and in contrast with that of a single multicomponent spin which turns around the energy 
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barrier. Interestingly, the Ising (hysteretic) behavior of the longitudinal component results from the interaction, and 
the presence of a transverse random field. 




FIG. 5. Energy profile of the two-spin toy model, for different values of the magnetic field H. At H — 0, the double well is 
symmetric. Larger values of H tilt the curve, eventually suppressing the left minimum. 



In the zero dimensional (single spin) or infinite dimensional (mean field) case, the system follows a single minimum 
and there is no hysteresis. When fluctuations are included, however, the above toy system shows that many minima 
appear, allowing for a hysteretic behavior. Thus, as in the scalar case, the action of Eq. (|l7|) has to be corrected in 
order to remove the unphysical minima. Following the procedure of Sec. IV. B, we replace s (x, t) by s (x, t) + <r (t), to 
obtain, if a is parallel to H, the corrected action 

S = S [s, s + a] 

= J dtd d xs\\ {-r}d t (a,, +a)+ KV 2 s\\ + H + c x («„ + a) + c 2 (s\\+a) 2 + s 2 ± (sy + a) H j 

/R 
dtdt'd d x — sy (x,i) sy (x,t') 



dtd d x s_l • | -rjd t s^ + K\7 2 s^ + CxS± + c 2 
+ J dtdt'd d x s± (x, t) ■ s_l (x, t') , 



(46) 



where the subscripts || and _L refer to the components parallel and perpendicular to the magnetic field. Expanding 
the polynomials, and absorbing the parameter a in thereby corrected coefficients, leads to the form 



dtd d 



-r]d t s\\ + XV 2 sy + A + Aisy + A 2 s\ + A 2 s 



2 1 A 3 sf + A 3 sis|| 



f R 

- j dtdt'd d x — sy (x, t) sy (x, t') 



dto fx s 



-r]d t s ± + KV 2 s ± + B 1 s ± + i? 2 S||Sj_ + -B 3 s^s ± + B 3 s 2 s ± 



/R 
dtdt'd d x — s_l (x, t) ■ s ± (x, t') . 



(47) 
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Clearly, only even powers of s± are present in the first line, while only odd powers occur in the third. A term 
of order p in s, trivially (i.e. ignoring the coarse graining step which couples different order terms) scales with 
b 2+d^ P = b 2 P -(p-i)d/2_ In dimensions d < 4, all terms are relevant, and a non-trivial fixed point at long length scales 
cannot in general be reached by tuning merely two quantities, namely the external magnetic field and the amount of 
randomness. Similarly for d = 5, in which terms up to p = 5 are relevant. In 6 — e dimensions, however, all terms 
with p > 3 are irrelevant under the RG, leading to an effective action of the form of Eq. (f47j), with any terms not 
displayed set to zero. Including the corrections due to coarse graining, to lowest non-trivial order (in the interaction 
and in e = 6 — d), the recursion relations for the remaining 10 vertices read 



A' 



b 3 ~^ 2 [A + (I + 2Ai) A 2 + (n 
b 2 ' 
b 2 



1) (/ + 2Bi) A 2 \ 
Ai + 3 (/ + 2Ax) A 3 + (n - 1) (J + 2Bi) A 3 ] 
B x + (n + 1) (J + 2Ai) B 3 + {I + 2B X ) B 3 ] 



A' n = 6 1+£ / 2 A 2 + 18A 2 A 3 + 2 (n - 1) A 2 B 3 
A 2 + 2A 2 A 3 + 2 (n + 1) A 2 B 3 + 
£ 2 + 4£ 2 B 3 + 2 (n + 1) S 2 B 3 + 
A 3 + ISA 2 + 2 (n - 1) A 3 B 3 ] 
A 3 + AA 2 3 + 6A 3 A 3 + 2(n + 1) A 3 B 3 
B 3 + 2 (n + 7) B| + 2A 3 B 3 ] 



= 

A' 3 = 
A' 3 = 
B< 3 = 
B' = 



b l+e/2 

If 
If 

V 



_3 

B 3 



ABl 



2 (n + 1) B 3j B 3 + 6A 3 B 3 



- 2 (n - 1) B 2 A 3 J 
2B 2 A 3 + 4A 2 A 3 ] 
2B 2 A 3 + 4A 2 /5 3 

+ 4A 3 B 3 ] 

+ 4A 3 B 3 ] 



4A 2 B 3 ] 



(48) 



where I = A 2 (b 2 - l) /26 2 In b, and a factor of R In 6/ (47r) is absorbed in a redefinition of A 2 , A 2 , i? 2 , A 3 , A 3 , B 3 , 
and i? 3 . (For the derivation of the corresponding relations for the scalar model, the reader is referred to Refs. 
the extension to the effective action of Eq. (E^) is straightforward.) 



2. Ising Criticality 

Clearly, a non-trivial fixed point for all 10 vertices cannot be reached in general if only two quantities are to be 
tuned. However, if R is appropriately tuned, A\ flows to its fixed (finite) value, while B\ grows indefinitely. Then, 
under the RG, the interactions in the perpendicular components become less and less important with respect to the 
quadratic term. After sufficient rcscalings, the theory becomes Gaussian in the perpendicular fields which can be 
integrated out, resulting in an effective action for S|| identical to the scalar action of Eq. (|39|). Thus, the critical point 
of our O (n) model is generically described by the same action as in Refs. ijl)-^, yielding identical recursion relations 
and exponents. This can be physically understood in the following way: if a configuration of several spins gives rise 
to many minima, a coarse-grained vector-spin system roughly looks like an Ising system, leading to scalar-like critical 
behavior. The latter occurs for a non-vanishing value of the magnetic field or the magnetization, at which the full 
rotational symmetry of the O (n) model is broken. 



3. Transverse criticality 

In the case considered above, the longitudinal field is massless and the transverse components massive. Alternatively, 
one should be able to reach another fixed point which incorporates the reverse situation. In Eq. (f47|), let us expand 
the longitudinal field about some parameter A (t), so chosen as to cancel the linear term Aq, i.e. we write 

S|| = s|| + A, (49) 

with A satisfying 

A + AiA + A 2 A 2 +A 3 A 3 = 0. (50) 

The theory is then Gaussian in s',, . Integrating out the longitudinal field, the action for the transverse field reduces 
to Eq. (|I^ ) with (§j_,Sj_) instead of (s, s), and Aq = A 2 = 0, i.e. 

S= I dtd d x s± ■ [(-ridt + KV 2 + r±) s± + B 3 s 2 ± s ± ] , (51) 
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where r± denotes the corrected mass (see below). Similarly to the scalar model, this action is critical for an appropriate 



value of the disorder, or equivalently of the "mass" . The recursion relations for r± and uj_ = R 
be read off from Eq. pq) by setting A l — A l = B 2 = B3 = 0, as 



In 6/ (4tt) 3 



r> x = b 2 {r ± + [{n-l) + 2}(I + 2r ± )u ± } 
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_L 



¥ { Ul _ + 2 [(n - 1) + 8] u 2 } 



B-i can 



(52) 



From these recursion relations we can obtain the static exponent v with which the correlation length £ diverges. By 

—i/ o ("- 1 )+ 2 C 

definition £ ~ (i? — i? c ) " ~ |r — r c | , which along with £' = £/b and (dr'/(9r) fixod point = 6 («-i)+s e = frs/^ yields 



(6 2/r ' |r — r c |) " = \r — r c \ " b , i.e. vy r — 1, and 

1 l(n-l) + 2 
v = - + -7 ; 

2 2(n-l) +8 



(53) 



Note that the action of Eq. ( |5l| ) is identical to the critical action for n — 1 (weakly coupled) scalar fields, with 
rotational symmetry. It follows immediately from this consideration that the recursion relations and exponents for 
the longitudinal (Ising-like) fixed point are identical to those calculated here, with n — 1 = 1. 
In the reduced action, the transverse components' bare mass becomes 

r ± = B 1 +B 2 X + B 3 X 2 . (54) 



The critical line r* x (u±) is in the lower half (r± < 0) of the (r±, u±_) plane. Since Aq — + ±00 as H ~ * ±00, Eq. (|5C 
implies that A also goes from — 00 to +00 as the field is increased (similarly to A3 < [Eq. (f42|)], we have B3 < 0). 
Thus, r± crosses the critical line for a given value of H, provided that B\ — 4B3 (B\ — r\ (u±)) > 0, i.e. unless B\ is 
too negative (for a Mexican hat potential, we expect B\ > 0). This may not be true if is a discontinuous function 
of H . But A (H) is clearly continuous for critical and higher disorders. (A very high disorder, though, may suppress 
this transverse criticality by renormalizing r± into large negative values.) Whether or not the fixed point controlling 
the transverse criticality is reachable experimentally depends on what region of the space of the theory's parameters 
is swept when the physical quantities at hand in the experiment are varied. Unlike many equilibrium problems in 
which the symmetry is broken by an infinitesimal field, criticality can occur here at large values of H, allowing for 
non- negligible higher order terms, such as H 2 s^ or H 2 s 2 s 2 , in the Hamiltonian. These terms appear in the theory as 

modified (strong) functional dependences of the coefficients Aj, A,-, Bi, Bi, on the magnetic field, which, depending 
on the trend, may favor a transverse instability. 

The transverse critical point corresponds to an infinite susceptibility at some H±, resulting in a spontaneous 
transverse magnetization. A similar phenomenon was noted for the case of a pure (thermal) system in Ref. [ [l6|| . 
There, however, the appearance of a transverse magnetization is due to a magnetic field which oscillates at high 
frequency. It is a purely dynamical effect, not observed in the 77 — > limit. In our case, the effect is due to the 
presence of a quenched randomness, and its non-equilibrium nature lies in the metastability of the minima occupied 
by the system during its history. This transverse instability, though, differs from the longitudinal criticality in that 
it occurs for a range of disorders, rather than at a specifically tuned amount of randomness. Its underlying physical 
mechanism may for example be illustrated by two beads, attached to each other by a spring, and flowing on the two 
sides of a Mexican hat's bump, as it is progressively tilted. At some point, it might become favorable for one of the 
beads to jump on the opposite side of the rim, and to continue its motion next to the other bead, corresponding to a 
transverse ordering. 

At H±, the transverse components choose one of the many minima, which breaks the rotational symmetry in the 
(71 — l)-dimensional transverse space. We therefore have to correct the action of Eq. (|47|) for H > H±, as we did 
for the longitudinal field. This however does not alter the analysis of the longitudinal criticality. By shifting sj_, we 
eliminate the transverse linear term and can then follow the same procedure as before. We used the fact that the 
correction to Eq. (E7|) vanishes for H — > ±00, and this clearly still holds here. 



4- O (n) criticality 

In the above, we examined two distinct fixed points, corresponding to an infinite longitudinal susceptibility and 
a transverse ordering, respectively. For a specific choice of the theory's parameters, the two should merge into a 
single, rotationally invariant (in the n-dimensional space of the fields) fixed point, at which all components become 
simultaneously critical. In addition to the magnetic field and disorder, how many quantities should we tune to reach 
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such an O (n) fixed point? While Eq. (K7T) is a suitable formulation of the model when the symmetry is broken in the 
longitudinal direction, the form of Eq. (|46|) is more appropriate close to a fully rotationally invariant theory. An RG 
transformation can be carried out for it, as was done for the action of Eq. fl47|). Since cr (t) is spatially uniform, only its 
q = mode is non-vanishing. The parameter cr, therefore, does not participate in the coarse graining transformation 
and its normalization is given by the rescalings of the coordinates and the fields, as 

o-'(t)=C 1 a(b 2 t)=bl- 2 o-(b 2 t). (55) 



The recursion relations for H, ci, and C2, on the other hand, are obviously given by Eq. (|4q), with Aq = H, 
Ai=B x = a, A 3 =A 3 =B 3 = B 3 = c 2 , and A 2 = A 2 = B 2 = 0. 

By tuning both the magnetic field H, and H , defined as the zero of a (a (H ) = 0), to zero, the action of Eq. ( |4^ ) 
reduces to the form of Eq. (|5l]), with n fields s rather than the n— 1 components of sj_, which is critical for a given 
value of the disorder. Consequently, the recursion relations and exponents are identical to those characterizing the 
transverse critical point, but with n replacing n — 1. Furthermore, since H = and the action is fully rotationally 
invariant, hysteresis is suppressed at the O (n) critical point; in particular, m (H = 0) = 0. Finally, by identifying the 
relevant parameter cr, we have established that only a single additional physical quantity needs to be tuned to reach 
the O (n) fixed point (provided that Hq crosses zero as the physical quantity in question is varied, at H = 0). 

The symmetric fixed point is characterized by three relevant parameters, namely r and H, with the usual exponents 
y r = 2 — e (n + 2) / (n + 8) and yn — 3 — e/2, and a with y a = d/2 — 2, to O (e). The parameter cr is a measure of 
the deviation from a symmetric (non-hysteretic) theory at H = 0, and, according to the above exponents, is much 
less relevant than the other symmetry breaking term H, or of the mass r, to 0(e). Ref. Q notes the fact that the 
"critical region" is unusually large in the scalar model, as signaled by power laws with surprisingly high cutoffs even 
a few percent to one hundred percent away from the critical disorder, and briefly discusses possible origins of this 
phenomenon. This observation, along with the weakness of cr's relevance, suggests in our case that even away from 
the O (n) fixed point, an O (n)-like behavior might be displayed at short ranges by the system, before crossing over 
to an Ising-like behavior at long enough time and length scales ||. 

As mentioned in the introduction, a vanishing magnetization does not a priori ensure a fully rotationally invariant 
system; higher moments could in general display anisotropy. In that case, a complete theory would break the rotational 
symmetry, leading to anisotropic exponents different from y r , yu, and y a above, which are the outcome of a restricted 
action that properly describes only quantities linear in the spin field. In general we may ask the following question, 
which applies to the three cases, Ising, transverse, and O (n). Would a complete theory yield the same exponents as 
our resticted theory? The latter correctly describes the magnetization m_ (t). That is, if one were able to calculate 
the path integral 

J VsVss(x,t)e s[g - s+a] , (56) 

one would obtain m_ as a function of H (t), R, and an additional tuning parameter (corresponding to Hq), and thus 
the associated exponents describing the critical singularity of the magnetization. It is easy to see that, as for the non- 
random equilibrium thermal Landau- Ginzburg model, the exponents y r , yu, and y a have a one-to-one correspondancc 
with the exponents that characterize the singular behavior of the magnetization, as well as with those associated with 
a number of other quantities, such as the correlation length. This immediately implies an affirmative answer to 
the above question; the requirement that the action should describe the physical magnetization is a strong enough 
constraint for the theory to correctly generate the singularity of, e.g., the correlation length. Evidently, this argument 
does not apply to the critical behavior of higher moments of the spin field distribution, and whether or not the latter 
are isotropic at criticality is an interesting open question. In fact, even the isotropy of y r at the O (n) critical point 
may seem surprising at first, since the history of the system apparently sets up an anisotropic context for the spins' 
motion. It is consistent, however, with the fact noted in Sec. IV. A, that different time slices of the action decouple 
under the RG in the r\ — > limit, in which the static exponents are calculated. 

As the reader may have noticed, our recursion relations and exponents in d — 6 — e are none other than the recursion 
relations and exponents for the pure (equilibrium) O (n) model in d = 4 — e. Indeed, dimensional reduction to two 
lower dimensions holds perturbativel y p| ]3T ]l^]i"5]| . This allows us to obtain, with no further calculational effort, the 
expansion for v to higher orders in e Jl8| . 

The dynamic exponent z, on the other hand, cannot be obtained from a time-independent model. How is the 
renormalization procedure modified when r\ ^ 0? In the coarse-graining process, vertices which couple different times 
are generated. In particular, a non-local quadratic term of the type 

t' 

dt Silvan (*)/(*'-*)» (57) 
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is generated in the action, where / is some function (containing free propagators, disorder, etc). Expanding su (t) 
about t' , this is rewritten as 



dt su (t') \ s|| (t') + 



dsn 
~8t 



(t-t') + ...\f{t'-t) 



(58) 



The first term contributes to the coarse-grained "mass" , the second results in a correction to z. In fact, the perturbative 
calculation we just briefly described is similar to that done by by Krey [[l9) for the time-dependent thermal random 
field model, where he calculates z to O (e 3 ), giving to lowest non-trivial order in e, 



n + 2 
(n + 8f 



e 2 . 



(59) 



(Here, n is the number of components that become massless at criticality.) 



V. CONCLUSION 



The main aim of the present paper was to extend earlier work on critical hysteresis to the case of an ordering of 
continuous symmetry. Our central result is that, generically, criticality is Ising-like, i.e. the critical exponents calcu- 
lated for a scalar field model still hold for a vector field. By tuning a single additional quantity, however, a fully 
symmetric fixed point can be reached, to which O (n)-likc exponents arc associated. Furthermore, the possibility of a 
spontaneous non-equilibrium transverse magnetization is unveiled and examined using a perturbative renormalization 
scheme. In addition, we have clarified several issues pertaining to the path integral formalism, in particular to the 
structure of time dependences, and to the problem of multiple solutions. These analyses are useful for our treatment 
of the problem, as well as for a clearer understanding of earlier works. 

An interesting question which remains open is that of the lower critical dimension. In the fully rotational case, 
naive dimensional reduction suggests that the lower critical dimension is 4. If this were the case, no vector criticality 
should be observed in three dimensions. Several scenarios are possible. For example, the hysteresis curve may display 
a jump for low disorder and be continuous for high disorder, but without the intermediate limit case of a continuous 
curve with a diverging slope at a given point. A more probable scenario is one in which the hysteresis curve is already 
smooth for an infinitesimal amount of disorder. 
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APPENDIX A: EXPANSION ABOUT MEAN FIELD THEORY 



In Refs. jl]-[|, the effective action is written as an expansion about mean field theory (MFT). In such a formulation, 
the linear term in the action (corresponding to Aq in Eq. (j47|)) is vanishing. The "masses", i.e. the coefficients r\\ and 
r± (corresponding to A\ and B\ in Eq. (|47|)) of the parallel and perpendicular quadratic terms are calculated as the 
elements of a particular response tensor within MFT ^^-^]. In this appendix, we calculate and r±, starting from 
the mean field equations of motion, and find that, generically, r« ^ r±. This appendix thus reaffirms that criticality 
is Ising-like in general, as found in Sec. IV, and makes contact with the earlier methodology of Refs. ffl-|3UTF 

MFT is defined by an infinite-range coupling, leading to the equations of motion 



rjd t s (x) = m (t) + H (t) + f (t) +h (x) + Cl s (x) +c 2 s (x) s (x), 



(Al) 



where m (t) is defined self-consistently by m(t) = s (x, t)| f _ , and f (t) is a test field. The "masses" are the static 
components of the tensor Ss/St — I, where / is the identity matrix p3|,p]-p| . We have 



Sm 



( 8m 1 (H+H 1 ,H 2 , 
1 5H X 
Sm 2 (H+H 1 ,H 2 , 
SHi 

V 



Sm 1 {H+H 1 ,H 2 , 
SH 2 

8yn 2 (H+H l ,H 2 , 
SH 2 



...\ 



Hi —H-2- 



( X|| (H) 

X± (H) 

xx (H) 

V • ; 



(A2) 
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(We consider 77 — > 0, so responses are non-vanishing only at equal times, and only the static part remains. Also, in 
MFT, the magnetization is always parallel to the magnetic field, leading to vanishing off-diagonal elements and to the 
relation x± = rn/H j|l.) In order to calculate the relevant tensor, let us split the field f between m and H, in such a 
way as to satisfy Eq. ( |A2|) for effective magnetization and magnetic field. That is, we write the equations of motion 
in the following fashion, with corrected m and H, and with no additional field, 



77<9 t S|| = \ m -+ 
for the longitudinal component, and 

r/d t s ± 



X|| 



1 + X|| ' 



,i + x± 

for the transverse components. Whence 



5s 
Si 



/II + H 



1 



/|| + h\\ + ciSj| + c 2 s s|| 



1 



1 + X-L 



fj_ + hj_ + C]S_L + C 2 S 2 Sj_ 































1 

i+x 



(A3) 



(A4) 



(A5) 



and 



/■ = , and r± = 



1 + X\\ 



l + X± 1 + m/H' 



(A6) 



Thus, rn (H) ^ r± (H) in general. 
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